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We report adaptive resolution molecular dynamics simulations of a flexible linear polymer in 
solution. The solvent, i.e., a liquid of tetrahedral molecules, is represented within a certain radius 
from the polymer's center of mass with a high level of detail, while a lower coarse-grained resolution 
is used for the more distant solvent. The high resolution sphere moves with the polymer and 
freely exchanges molecules with the low resolution region through a transition regime. The solvent 
molecules change their resolution and number of degrees of freedom on-the-ny. We show that 
our approach correctly reproduces the static and dynamic properties of the polymer chain and 
surrounding solvent. 
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I. INTRODUCTION 

The structure of polymers in solution is determined 
by the solvent-polymer interaction. In the case of a 
nonpolar polymer in an nonpolar solvent, one typi- 
cally distinguishes three "types" of solvent, good, or 
marginal, and poor. In the case of a good solvent the 
solvent-solvent, solvent-polymer, and polymer-polymer 
interactions effectively result in a situation, where the 
chain monomers are preferably surrounded by solvent 
molecules. As a result the chains are extended and the 
size scales as (R 2 oc N 2v ) with N being the number of 
monomers and v = 0.6 in three dimensions. For poor 
solvent one observes just the opposite and the chains col- 
lapse into a dense globule, (R 2 oc TV 2 / 3 ). The 9 regime 
is where these two effects compensate and the chains be- 
haves to a first approximation as a random walk, i.e., 
(R 2 cx N). In the limit of N — > oo the 0-point is a 
tricritical point in the phase diagram. As long as the sol- 
vent does not induce special local correlations beyond an 
unspecific attraction/repulsion and one is not studying 
dynamical properties the collapse of polymers is usually 
studied with an implicit solvent. The complicated local 
interactions are accounted for by an effective interaction 
between the chain beads. Studies of that kind have a long 
tradition in polymer science and the behavior is now well 
understood. Beyond that there are however many situa- 
tions, where it becomes difficult or even questionable, to 
ignore the local structure of the solvent. Solvent can play 
an important role in the functional properties of macro- 
molecules. For example, dehydration studies of proteins 
solvatcd in water demonstrated that at least a monolayer 
of water is needed for full protein functionality!!]. The 
influence of a macromolecule on the structure and dy- 
namics of the surrounding solvent is also an important 
issue. Therefore, a detailed study of interactions of a 
macromolecule with a solvent beyond effective coupling 
parameters is quite often required for an understanding 
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of the macromolecule's structure, dynamics, and func- 
tion. To determine the interactions of a solvent with a 
macromolecular solute chemistry specific interactions on 
the atomic level of detail have to be considered. However, 
the resulting solvating phenomena manifest themselves at 
mesoscopic and macroscopic scales [2[ and in the overall 
structure of the chains. Due to large number of degrees 
of freedom (DOFs) such systems are difficult to tackle us- 
ing all-atom computer simulations Q . Moreover, the vast 
majority of the simulation time is typically spent treat- 
ing the solvent and not the polymer or protein. A step to 
bridge the gap between the time and length scales acces- 
sible to simulations that still retain an atomistic level of 
detail and the solvating phenomena on longer time and 
larger length scales, is given by hybrid multiscalc schemes 
that concurrently couple different physical descriptions of 
the system (see, e.g., Refs. 0, H, ©, 0, fi ©, M ED] ) ■ 

Recently, we have proposed an adaptive resolution 
molecular dynamics (MD) scheme (AdResS) that concur- 
rently couples the atomistic and mesoscopic length scales 
of a generic solvent [H, EH, G3 ■ In the first application we 
studied a liquid of tetrahedral molecules where an atom- 
istic region was separated from the mesoscopic one by a 
flat or a spherical boundary. The two regimes with dif- 
ferent resolutions freely exchanged molecules while main- 
taining the thermodynamical equilibrium in the system. 
The spatial regions of different resolutions, however, re- 
mained constant during the course of the simulations. 
More recently this approach was extended to the study 
of water [l5j . In the present paper, we generalize our 
approach to the study of a polymer chain in solution. 
The chain is surrounded by solvent with " atomistic" res- 
olution. When the chain moves around, the sphere of 
atomistically resolved solvent molecules moves together 
with the center of mass of the chain. In this way the 
chain is free to move around, although the explicit res- 
olution sphere is much smaller than the overall simula- 
tion volume. This enables us to efficiently treat solvation 
phenomena, because only the solvent in the vicinity of 
a macromolecule is represented with a sufficiently high 
level of detail to take the specific interactions between 
the solvent and the solute into account. Solvent farther 
away from the solute, where the high resolution is no 
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longer required, is represented on a more coarse-grained 
level. In this work, a macromolecule is represented by 
a generic flexible polymer chain embedded in a sol- 
vent of tetrahedral molecules introduced in Refs. [I5,ll3l|. 
This study represents a first methodological step towards 
adaptive resolution MD simulations of systems of biolog- 
ical relevance, e.g., a protein in water. 

The paper is organized as follows: In section II the dual 
scale model of a polymer chain in a liquid is presented. 
The hybrid numerical scheme and computational details 
are given in section III. The results and discussion are re- 
ported in section IV, followed by a summary and outlook 
in section V. 



II. MULTISCALE MODEL 

We study a single generic bead spring polymer solvated 
in a molecular liquid as illustrated in figure [TJ Solvent 





FIG. 1: (Color online) A schematic plot of a solvated generic 
bead-spring polymer. The solvent is modeled on different lev- 
els of detail: solvent molecules within a certain radius from 
the polymer's center of mass are represented with a high 
(atomistic) resolution while a lower mesoscopic resolution is 
used for the more distant solvent. The high resolution sphere 
moves with the polymer's center of mass. The polymer beads 
are represented smaller than the solvent molecules for presen- 
tation convenience, for details see text. 

molecules within a distance ro from the polymer's center 
of mass are modeled with all 'atomistic' details to prop- 
erly describe the specific polymer-solvent interactions. 
For the description of the solvent farther away, where 
the high resolution is not required, we use a lower reso- 
lution. The solvent molecules then, depending on their 
distance to the polymer's center of mass, automatically 
adapt their resolution on-the-fly. 

The model solvent is a liquid of n tetrahedral molecules 
as introduced in Refs. [l2l. Il3||. The solvent molecules in 
the high resolution regime are composed of four equal 
atoms with mass tuq. Their size a is fixed via the repul- 



sive Weeks-Chandler- Andersen potential 
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with the cutoff at 2 1 / 6 c a and e are the standard 
Lennard Jones units for lengths and energy respectively. 
fiaj/3 is the distance between the atom i of the molecule a 
and the atom j of the molecule /3. The neighboring atoms 
in a given molecule a are connected by finite extensible 
nonlinear elastic (FENE) bonds 
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with divergence length R = 1.5(7 and stiffness k — 
30e/cr 2 , so that the average bond length is approximately 
0.97cr for ksT — e, where T is the temperature of the 
system and ks is Boltzmann's constant. For the coarse- 
grained solvent model in the low resolution regime we 
use one-site spherical molecules interacting via an effec- 
tive pair potential [3], which was derived such that the 
statistical properties, i.e., the center of mass radial dis- 
tribution function and pressure, of the high resolution 
liquid are accurately reproduced. This is also needed for 
the present study, since the motion of the high resolution 
sphere should not be linked to strong rearrangements in 
the liquid. The high and low resolution freely exchange 
molecules through a transition regime containing hybrid 
molecules (see figure [1} , where the molecules with no ex- 
tra equilibration adapt their resolution and change the 
number of DOFs accordingly p^. [l3l [iH . 

The polymer is modeled as a standard bead-spring 
polymer chain [l6j. It contains N monomers, which rep- 
resent chemical repeat units, usually comprising several 
atoms. The interactions between monomers (beads) are 
defined using Eqs. fT} and @ with the rescaled val- 
ues ob — 1.8(7, Ro B = RqObIo = 1-5(7.b, and ks = 
ka 2 /a B = 30e/a B , such that the size of the polymer bead 
<jb is approximately the same as the size of the solvent 
molecule [13 ■ The average bond length between beads 
is rescaled accordingly. The bead mass is also increased 
rriB = 5mo to make them behave more like Brownian 
particles. Standard Lorentz-Berthclot mixing rules [l7| 
are used for the interaction between monomers and the 
'atoms' of the solvent molecules. 



III. NUMERICAL SCHEME AND 
COMPUTATIONAL DETAILS 

To smoothly couple the regimes of high and low level 
of detail of the description of the solvent molecules, we 
apply the recently introduced AdResS scheme [l2j. There 
the molecules can freely move between the regimes, they 



3 



are in equilibrium with each other with no barrier in be- 
tween. The transition is governed by a weighting function 
w(r) e [0, 1] that interpolates the molecular interaction 
forces between the two regimes, and assigns the identity 
of the solvent molecule. We resort here to the weighting 
function defined in Rcf. Hal: 



w(r) 
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[Ur-r )] 



r > r > 
r > ro + d 
r + d > r > r 



(3) 



where ro is the radius of the high resolution region and 
d the interface region width, cf. [1] The radius ro must 
be chosen sufficiently large so that the whole polymer 
always stays within the high resolution solvent regime. 
w(r) is defined in such a way that w — 1 corresponds 
to the high resolution, w = to the low resolution, and 
values < w < 1 to the transition regime, respectively. 
This leads to intermolecular force acting between centers 
of mass of solvent molecules a and 0: 
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F a p is the total intermolecular force acting between cen- 
ters of mass of the solvent molecules a and (3. F^j 
is the sum of all pair 'atom' interactions between ex- 
plicit tetrahedral 'atoms' of the solvent molecule a and 
explicit tetrahedral 'atoms' of the solvent molecule j3, 
F^o is the effective pair force between the two solvent 
molecules, and R a , R^, and R are the centers of mass of 
the molecules a, (3 and the polymer, respectively. Note 
that one has to interpolate the forces and not the in- 
teraction potentials in Eq. if the Newton's Third 
Law is to be satisfied 14]. To suppress the unphysi- 
cal density and pressure fluctuations emerging as arti- 
facts of the scheme given in Eq. ((4]) within the transi- 
tion zone we employ an interface pressure correction [l3j. 
The latter involves a reparametrization of the effective 
potential in the system composed of exclusively hybrid 
molecules with w = 1/2. Each time a solvent molecule 
crosses a boundary between the different regimes it gains 
or looses on-the-fly (depending on whether it leaves or 
enters the coarse-grained region) its equilibrated rota- 
tional and vibrational DOFs while retaining its linear 
momentum fLU [Lsl [T^. |20T] . This change in resolution re- 
quires to supply or remove "latent heat" and thus must 
be employed together with a thermostat that couples lo- 
cally to the particle motion [l2|, EH ■ This is achieved by 
coupling the particle motion to the Dissipative Particle 
Dynamics (DPD) thermostat (2lj. This bears the addi- 
tional advantage of preserving momentum conservation 
and correct reproduction of hydrodynamics in our nVT 
MD simulations. Because of the freely moving polymer 
chain and solvent molecules, the above scheme requires 
the center of the high resolution sphere to move with the 
polymer but slowly compared to the surrounding solvent 
molecules, so that they at the boundary between differ- 
ent regimes have enough time to adapt to the resolution 



change. The validity condition for our approach thus re- 
quires D polymer *^ ^solvent ; where Dpolymer and D solvent 

the corresponding diffusion constants. This condition is 
trivially fulfilled in polymeric solutions and thus also in 
our simulations (see the next section). 

We conducted all MD simulations using the ESPResSo 
package [22j. We integrated Newtons equations of mo- 
tion by a standard velocity Verlet algorithm with a time 
step At = 0.005r and coupled the motion of the parti- 
cles to a DPD theromstat[21J with the temperature set 
to T = e/k B - The DPD friction constant C = 0.5T" 1 , 
where r = (e/niQij 2 )^ 1 ^ 2 , and the DPD cutoff radius was 
set equal to the cuttoff radius of the effective pair interac- 
tion between solvent molecules, i.e., 3.5cr[l3j]. The width 
of the transition regime is 2.5c[l~3| . Periodic boundary 
conditions and the minimum image convention [17[ were 
employed. After equilibration, trajectories of 5000r were 
obtained, with configurations stored every 5t. These pro- 
duction runs were performed with a 10 9 e/a force capping 
to prevent possible force singularities that could emerge 
due to overlaps with the neighboring molecules when a 
given molecule enters the transition layer from the coarse- 
grained side (l2| . The temperature was calculated using 
the fractional analog of the equipartition theorem: 
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where (K a ) is the average kinetic energy per fractional 
quadratic DOF with the weight w(r) = a[l4(. Via 
Eq. (|5|) the temperature is also rigorously defined in 
the transition regime in which the vibrational and rota- 
tional DOFs are partially 'switched on/off'. The molec- 
ular number density of the solvent is p = 0.175/<t 3 , 
which corresponds to a typical high density Lennard- 
Jones liquid[13l|. We considered three different sys- 
tem sizes with corresponding cubic box sizes: L = 
25.0a, 30.6cr, 34.2<r. The reduced Lennard- Jones units (l7j 
are used in the remainder of the paper. 



IV. RESULTS AND DISCUSSION 

To validate the AdResS approach for the present poly- 
mer solvent system, we carried out the analysis of the 
structural and dynamic properties of a polymer chain 
in the hybrid multiscale solvent compared to the corre- 
sponding fully explicit system where all solvent molecules 
are modeled with a high level of detail, i.e., as a tetrahe- 
dral molecules. 



A. Statics of the Polymer Chain and Solvent 

First, we focus on the explicit (ex) systems where the 
solvent is modeled with the high resolution all over the 
simulation box. These results are considered as the refer- 
ence to check how well AdResS produces the same physics 
as the all-atom MD simulation. The reference average 
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thermodynamic properties of the corresponding ex sys- 
tems (polymer+explicitely resolved solvent) are listed in 
table III 



N 


10 


20 


30 


L 


25.0 


30.6 


34.2 


<p> 


2.01 ±0.04 


2.01 ±0.03 


2.02 ±0.01 


<T> 


1.0 ±0.01 


1.0 ±0.01 


1.0 ±0.01 


^ Tp ly mer ^> 


1.0 ±0.5 


1.0 ±0.3 


1.0 ±0.2 


Tsolvent 


1.0 ±0.01 


1.0 ±0.01 


1.0 ±0.01 



TABLE I: Thermodynamic properties of the fully explicit sys- 
tems (ex) of different chain lengths TV and box sizes L: aver- 
age total pressure (p), average total temperature of the sys- 
tem (T), average temperature of the polymer (T po i ymer ) , and 
average temperature of the solvent (T so i vent ). 

The static properties of the solvent are characterized 
by the solvent radial center of mass distribution (RDF) 
function depicted in figures [2] (a). This distribution func- 
tion is within the thickness of the lines the same for all 
systems studied, including the hybrid ones. This is to 
be expected from our previous studies and the fact, that 
the polymer fraction of volume is very small compared 
to that of the solvent. 

The statistical properties of polymers are conveniently 
described by a number of quantities, namely the radius 
of gyration 

<^> = ^E<( r «- R ) 2 >> (6) 

i 

where is the position vector of the ith monomer and 
R = /V -1 J^i r i is t ne polymer's center of mass, the end- 
to-end distance 

(R 2 E ) = ((r N -r 1 ) 2 ), (7) 
and the hydrodynamic radius 




where r i} = |r< — r j | [23l| . 

(R%) and (R%) scale as 

(R 2 g )k(R%)<xN 2 * (9) 

with the number of monomers N where v = 0.5 in 9 sol- 
vent and v « 0.588 in good solvent conditions [24l l25l l2r| 
with rather small finite size corrections, while the hydro- 
dynamic radius is known to show significant deviations 
from asymptotic behavior up to very long chains [27j. 
The single-chain static structure factor S(q) 
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FIG. 2: (a) The solvent center of mass RDF for three dif- 
ferent polymer lengths (N=10, 20,30), (b) the static structure 
factor of the polymer in the Kratky representation, and (c) 
the solvent density around the center of mass of the chains, 
which illustrates the so called correlation hole. 

probes the self similar structure within the scaling regime 
and thus provides an accurate way to determine v. S(q) 
scales as 

S(q) OC q- 1 '" -> q 2 S(q) oc q 2 ' 1 ^ (11) 

in the regime Rq 1 -C q <C b^ 1 , where b is the typical 
bond length. By fitting a power law to the computed 
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q 2 S(q) plotted in figure [2] (b) we obtained the values for 
v reported in table |TTJ Table |TT] summarizes the values of 
all quantities defined above, which characterize the static 
properties of the polymer chain. The calculations were 
performed for N = 10, 20, 30. 



N 
L 



Rg = {RlY^ 

Re = \R%) 

Rh = (RhT 1 

V 

2/z 

r = R 2 e/(6D) 



10 


20 


30 


25.0 


30.6 


34.2 


4.2 ±0.8 


5.7 ± 1.0 


8.1 ± 1.4 


2.7 ±0.5 


3.8 ±0.6 


5.0 ±0.8 


6.7 ±2.0 


8.6 ±3 


12 ±3 


3.3 ±0.3 


4.0 ±0.3 


4.7 ±0.4 


0.63 


0.54 


0.57 


0.59 


0.71 


0.67 


152 


481 


1390 



TABLE II: Summary of some polymer (embedded in the 
explicitely resolved ex solvent) properties: number of poly- 
mer beads N, size of the simulation box L, average maximal 
distance of a monomer from the polymer's center of mass 
{Ar m ax), radius of gyration Rg, end-to-end distance Re, hy- 
drodynamic radius Rh, the static exponent v, the exponent 
2/z, where z is the dynamic exponent, and the longest relax- 
ation time t calculated using data from table PVTl The error 
bar for the exponents v and 2/z is roughly 10%. 

Another property, which directly reveals the fractal 
structure of the chains is the correlation hole, which is 
shown in figure[2](c). It directly shows, to which distance 
from the center of mass of the chains, the solvent density 
is perturbed by the chain beads. For the later applica- 
tion of the hybrid scheme it is important to define the 
explicit solvent regime large enough in order to cover the 
correlation hole completely. 

The values of v actually differ slightly from the asymp- 
totical value for the good solvent due to the finite chain 
lengths. Nevertheless, the agreement improves with the 
increasing N, as expected. 

Let us now turn our attention to the hybrid solvent 
studied by MD simulation using AdResS. 

To assure that the polymer is surrounded only by 
the explicitely resolved molecules, we determine first the 
maximal monomer distance from the polymer's center of 
mass, Ar ma x, as shown in figure [3] for all chain lengths 
studied. As shown Ar max always stays within the high 
resolution regime. 

Because for N = 30 Ar max gets rather close to r , 
we checked the static polymer properties for that case 
again. In figure 3] we compare the all explicit simulation 
to the two hybrid simulation schemes (with (ex-cgi C ) and 
without (ex-cg) the pressure correction [13J in the transi- 
tion regime) for the chain form factor and the correlation 
hole. The agreement is excellent, showing that the pro- 
posed scheme should at least be capable of properly re- 
producing the conformational statistics of the embedded 
polymer in solution. 

This is first checked by comparing the average thermo- 
dynamic properties as given in table Hill for the hybrid ex- 
cg and ex-cgi C systems (polymer±hybrid solvent). While 
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FIG. 3: (Color online) Time evolution of the maximal 
monomer distance from the polymer's center of mass, Ar m ax, 
for polymers with N = 10 beads and the radius of the high 
resolution regime ro = 7.0 (red line), N = 20 and ro = 11.0 
(green line), and N = 30 and rg — 12.0 (blue line). 



10 



o.i 



1 1 


1 1 1 1 1 1 1 1 1 / 1 1 1 1 1 1 




ex N =311 




ex - cg N=ia ■ 




ex-cg icN=30 

1 



0.1 



10 



(a) 



1.4 

I. 2 

1 - 
0.8 

II. 
0.4 h 
0.2 



1 r 

exw=30 
ex - c5jv =30 
ex — cgi CN ^ 3Q 




12 14 



10 



(b) 



FIG. 4: (a) The the static structure factor of the polymer with 
N = 30 in the Kratky representation for all three cases stud- 
ied: the fully explicite, the AdResS scheme with and without 
the interface pressure correction, (b) The correlation hole for 
the same systems as in (a). 



the temperatures are identical the pressure correction in 
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the interface layer reduces the pressure slightly, so that 
the hybrid system now also there agrees quite well with 
the all explicit simulation. The agreement with the ref- 



N 


10 


20 


30 


L 


25.0 


30.6 


34.2 


^ P ^ex — cg 


2.03 ±0.02 


2.04 ±0.01 


2.04 ±0.03 


*v P ^ &x — cgi c 


2.01 ±0.01 


2.01 ±0.01 


2.01 ±0.01 


< T > 


1.0 ±0.02 


1.0 ±0.01 


1.0 ±0.01 


^ '^■polyme.r ^ > 


1.0 ±0.5 


1.0 ±0.3 


1.0 ±0.2 


^ solvent ^> 


1.0 ±0.02 


1.0 ±0.01 


1.0 ±0.01 



TABLE III: Thermodynamic properties of systems with the 
polymer solvated in the hybrid ex-cg solvent and the hy- 
brid ex-cgi c solvent: average total pressure (p), average to- 
tal temperature of the system (T), average temperature of 
the polymer (T po i ymer ), and average temperature of the sol- 
vent {T so i ven t). For the temperatures the results cannot be 
distinguished. 

erence values from table Q] is very good. This is in line 
with the general static properties of the polymers, which 
are given in llVl andlVl and compare very well to the data 
from tabic |TTJ 



N 


10 


20 


30 


L 


25.0 


30.6 


34.2 


ro 


7.0 


11.0 


12.0 


(^fmax ) 


4.2 ±0.8 


6.6 ± 1.2 


7.7 ± 1.3 


Rg = (R 2 a) 1/2 
Re = <i?i) 1/2 

r h = (R-yy 1 


2.7 ±0.4 


4.0 ±0.6 


4.6 ±0.7 


6.7 ±2.0 


10.4 ±2.8 


10.8 ± 2.5 


3.3 ±0.3 


4.0 ±0.3 


4.5 ±0.4 


V 


0.63 


0.58 


0.54 


2/z 


0.56 


0.69 


0.62 


r = R 2 G /(6D) 


122 


381 


882 



TABLE IV: Summary of some polymer (embedded in the hy- 
brid ex-cg solvent) properties: number of polymer beads N, 
size of the simulation box L, radius of the high resolution 
regime ro, average maximal distance of a monomer from the 
polymer's center of mass (Ar ma x), radius of gyration Rg, 
end-to-end distance Re, hydrodynamic radius Rh, the static 
exponent v, the exponent 2/z, where z is the dynamic expo- 
nent, and the longest relaxation time r calculated using data 
from table IVll The error bar for the exponents v and 2/z is 
roughly 10%. 

From the presented results we can conclude that 
AdResS faithfully reproduces the reference statics ob- 
tained from the simulations with a polymer embedded 
in the explicitely resolved solvent. 



B. Dynamics of the Polymer Chain and Solvent 

While the conformational properties of the polymer in 
solution are well understood and properly described by 
the adaptive resolution approach, the situation for the 
dynamics is much less clear. By changing the degrees 



N 


10 


20 


30 


r 

L 


25.0 


30.6 


34.2 


ro 


1 .0 


11. U 


ion 
12. U 




4.U ± U.8 


0.9 ± 1.2 


o a _l_ 1 k 
O.D ± 1.0 


Rg = (Ra) 1 


2.7 ±0.5 


3.9 ±0.7 


5.2 ±0.8 


D , 2 \l/2 

Re = (R E ) 


6.6 ± 2.1 


9.4 ± 3.0 


13.3 ± 3.3 


Rh = (R-yy 1 


3.2 ±0.3 


4.0 ±0.4 


4.8 ±0.4 


V 


0.59 


0.55 


0.57 


2/z 


0.69 


0.71 


0.77 


r = R 2 G /(6D) 


152 


362 


1127 



TABLE V: Same data as in table HVl but now for the hybrid 
ex-cgi c , where interface pressure correction is applied. 

of freedom not only the structure but also the dynam- 
ical properties are altered, however in a way which is 
less understood. It is also not a priori clear, whether an 
approach, which produces a precise coarse graining for 
structural properties, does this for dynamical properties 
as well. In a recent study of small additive molecules to a 
polymer melt it was shown, that while the length scaling 
is identical, the time scaling can be different [28| • In the 
present situation the influence of the transition regime 
poses additional difficulties. 

In order to determine the dynamical properties of the 
solvent and solute we calculated the respective diffusion 
coefficients. The diffusion coefficient of a species is com- 
puted from the center of mass displacements using the 
Einstein relation 

D= l lim (lR^)-R,(0)H = l 1 , m (Ag ! ) i 

6 t— >oo t 6 t^oo t 

where (t) is the center-of-mass position of the molecule 
i (which can be either a solvent or a solute molecule) 
at time t and averaging is performed over all choices of 
time origin and, in the case of solvent, over all solvent 
molecules. 

Figure [5] shows this for the solvent molecules' centers 
of mass as a function of time for the different systems 
indicated. All the curves in figure [51 except the one for 
the coarse-grained solvent coincide. Thus the effect of the 
polymer on the diffusivity of the solvent molecules is neg- 
ligible. In other words, the dilution is strong enough that 
the polymer effect on the solvent dynamics is very small. 
The coarse grained solvent molecules however move faster 
than the explicit ones. This is a consequence of the re- 
duced number of DOFs causing a time scale difference in 
the dynamics of the coarse-grained svstem[T^.[T3l|. While 
this can be very advantageous in some cases, one can 
also adjust D by an increased background friction in the 
DPD thermostat [U, HI- The diffusion coefficient D of 
the solvent was obtained by fitting a linear function to 
the curves depicted in figure[5j(a) and the obtained values 
are D bu i k(s:c = 0.036 and D bu i kcg = 0.057 for the explicit 
and coarse-grained solvent, respectively. The question to 
ask here is, to what extend does this have any influence 
on the dynamics of the embedded polymer. 
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FIG. 5: Log-log plot of the time dependence of the mean 
square displacement of the solvent molecule's center of mass 
in the. time interval [0,5000]: explicitely resolved (bulk e:r ) 
and coarse-grained (bulk C9 ) solvents without solvated poly- 
mer and the explicitely resolved solvent for the systems with 
three different lengths of the solvated polymer (N=10,20,30). 

Within the Zimm model[25j for polymer chain dynam- 
ics, which is known to describe the scaling of the dynam- 
ics in dilute solutions of polymers rather well and which 
takes into account the hydrodynamic interactions, the 
polymer diffusion coefficient scales as 

D oc N~" oc R H X oc R g \ (13) 

The longest relaxation time r = Rq/(6D), i.e., the Zimm 
time Tz oc R^ — R Z G , is the time the chain needs to move 
its own size, z = 3 is the dynamic exponent. Note that 
the motion of inner monomers within the appropriate 
scaling regime should be independent of N. For the mean 
square displacements of the monomers a scaling analysis 
immediately yields for the mean square displacement of 
a monomer i, 

{Ar 2 ) = (( Ti (t) - r 4 (0)) 2 ) oc t 2 l z = t 2 l\ (14) 

for distances significantly larger than the bond length and 
smaller than (i? 2 ), i.e., times smaller than tz- For the 
center of mass of the chains a diffusive behavior for the 
mean square displacement (AR 2 (t)) is always observed. 
Although the chains are relatively short, at least for N — 
30 one expects a behavior relatively close to the above 
mentioned idealized scheme [27j . Figure [6] shows (Ar 2 (t)) 
and (AR 2 (t)) for polymer chains with N = 10,20,30 
embedded in the different solvent scenarios studied. 

In all cases the observed exponents for (Ar 2 (t)) and 
(AR 2 (t)) are close to 0.7 ±0.05 and 1 ±0.05 respectively. 
Also the amplitudes of the displacements of the inner 
beads are almost the same for all chain lengths stud- 
ied. This is in good agreement with earlier studies on 
different generic polymer models in a explicit solvent as 
well as studies of chains in a hybrid lattice Boltzmann 
solvent HH2I. 

This is to be expected since we preserve 
the hydrodynamic interactions by employing the DPD 




(c) 

FIG. 6: Log-log plot of the time dependence of the mean 
square displacement of a single monomer (considered are only 
monomers near the chain's center of mass) for polymers and 
their centers of mass with N = 10, 20, 30 solvated in the ex 
solvent (a), the hybrid ex-cg solvent (b) and the ex-cgi C sol- 
vent (c) as indicated. 



thermostat in our simulations. This suggests that our 
hybrid scheme is also applicable to study dynamic prop- 
erties of a polymer in a solution. Small deviations how- 
ever occur in the diffusion constant itself. The diffusion 
constants for the polymers with N = 10, 20, 30 in the hy- 
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brid ex-cg and ex-cgi C solvents were obtained by fitting 
the straight curve to the polymer's center of mass mean 
square displacement presented in figure [6] The fit yields 
data listed in table IVIl The corresponding static and 
dynamic exponents and the longest relaxation times are 
given in tables IIVI and [V] 



N 


10 


20 


30 


D(ex) 
D(ex-cg) 
D(ex-cg ic ) 


0.008 
0.009 
0.0085 


0.005 
0.006 
0.006 


0.003 
0.0045 
0.0035 



TABLE VI: Diffusion constant of the polymer chain embed- 
ded in three different solvents: explicitely resolved ex, hy- 
brid ex-cg, and hybrid ex-cg ic . Though the statistics of the 
data is rather poor we can estimate the error bar roughly to 
10 — 15%. For comparison, the diffusion coefficients of the 
explicit and coarse-grained solvents are Db u ik C3: = 0.036 and 
Dbuik cg ~ 0.057, respectively with an effect of the polymers 
too small to determine here. Hence D po i ymer <C D so i vent - 

While the ratio of the diffusion constants for different 
chain lengths roughly follow the expected scaling, even 
though it cannot hold precisely due to the different box 
sizes, we here observe a tendency to a weakly accelerated 
diffusion in the hybrid regime. This is most evident for 
the hybrid ex-cg case. Two different aspects might play a 
role here. First the viscosity in the coarse grained outer 
regime is smaller, which must have an effect on the diffu- 
sion. Second, the small pressure and density fluctuations 
in the transition regime might contribute to the effect as 
well. 

Although this is only a very first and incomplete test, it 
shows that within the AdResS scheme essential aspects of 
the dynamical properties of the embedded polymer chain 
are reasonably well reproduced. 



V. SUMMARY AND OUTLOOK 

In this paper we presented a hybrid multiscale MD 
simulation of a generic macromolecule in a solvent us- 
ing the recently proposed AdResS method. The solvent 
surrounding the macromolecule is represented with a suf- 
ficiently high level of detail so that the specific interac- 
tions between the solvent and the solute are correctly 
taken into account. The solvent farther away from the 
macromolecule, where the high resolution is not needed, 
is represented on a coarse-grained level. The high and 
low resolution regimes freely exchange solvent molecules, 
which change their resolution accordingly. To correctly 
simulate momentum transport through the solvent, we 
use the DPD thermostat. The simulation results show 
that AdResS accurately reproduces the thermodynamic 
and structural properties of the system. The presented 
methodology is an extension of AdResS to simulations 
of a solvation cavity, and represents a first step towards 
the treatment of more realistic systems such as biomacro- 
molecules embedded in water. Work along these lines is 
already under way [To]]. 
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